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Abstract 

^~^ [ We reexamine, within the nonrelativistic constituent quark model (NRCQM), a recent claim that the current data 

'— ' on the strangeness form factors indicates that the uudss component in the proton is such that the uuds subsystem has 

Z^ ^the mixed spatial symmetry [31] ^y and flavor spin symmetry [4]i?5[22]f^[22]5, with s in S state (configuration I). We find 

■this claim to be invalid if corrected expressions for the contributions of the transition current to G^ and G|; are used. 

rr^We show that, instead, it is the lowest-lying uudss configuration with uuds subsystem of completely symmetric spatial 

ry^ -symmetry [4]^: and flavor spin symmetry [4]i?s[22]i?[22]s, with s in P state (configuration If), which could account for 

|the empirical signs of all form factors GJ;, G\j, and G\. Further, we find that removing the ccntcr-of-mass motion of the 

, clusters will considerably enhance the contributions of the transition current. We also demonstrate that it is possible 

|to give a reasonable description of the existing form factors data with a tiny probability Pss ~ 0.025% for the uudss 

I""! component. We further see that with a small admixture of configuration I , the agreement of our prediction with the 

'"pi data for G^ at low-g^ region can be markedly improved. We find that without removing CM motion, Pss would be 

overestimated by about a factor of four in the case when transition current dominates. We also explore the consequence 

'of a recent estimate reached from analyzing the existing data on d — u , s + s , and u + d — s — s, that Pss lies between 

•2.4 — 2.9%. ft would lead to a large size for the five-quark system and a small bump in both G|; 






f]G%j and G|; in the 
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region of q < 0.1 GcV within the considered model. 
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I The first indications of possible existence of strangeness 
content in the proton came from deep-inelastic muon scat- 
jtering, elastic neutrino-proton scattering, and analyses of 
■ttA^ cr-term [1|. Many other observables were later sug- 
gested, including excess 4> production in pp annihilation 
■[2], double polarizations in photo- and clcctroproduction 
'of (j) meson [3|, and asymmetry in scattering of longitu- 
dinally polarized electrons from polarized targets. While 
|the measurement of double polarization in (jj photoproduc- 
ition is being pursued with the development of polarized 
HD target by LEPS at SPring-8 [J], four vigorous experi- 
mental programs SAMPLE @, HAPPEx [g, A4 0, and 
GO [8| have already been undertaken to measure parity- 
violating asymmetry of polarized electron-proton scatter- 
ing in order to extract proton strangeness electromagnetic 
form factors. 

On the theoretical side, lattice QCD remains the only 
theoretical method which could provide a reliable deter- 
mination of the strangeness form factors from the first 
principle. For example, a recent low-mass quenched lat- 
tice QCD simulation, with the use of chiral extrapolation 
technique and the assumption of charge symmetry, gives 
^s = (-0.046 ± O.G19)Aijv S and G%{Q^ = 0.1 GeV^) = 



-0.009 ± 0.005 ± 0.003 ± 0.027 [l^. More recent LQCD 
efforts can be found in Ref. [IjJ. Nevertheless, a study of 
this intriguing question within hadron models could still 
provide invaluable insight concerning the underlying quark 
structure, e.g., how the strange quarks are arranged inside 
proton. 



Recently, Zou and Riska [12| considered the possible 
low-lying configurations of the uudss component of the 
proton within a constituent quark model. It was concluded 
that the empirical indications of a positive strangeness 
magnetic moment of the proton [5|, la, III suggest that the 
dominant ss configuration in the proton would have the 
s in the ground state and uuds system in the P state. 
It would lead to an interesting implication that the qqqqq 
components in the proton would mainly be in colored quark 
cluster configurations rather than in "meson cloud" con- 
figurations as commonly perceived. The calculation of 
Ref. [12| was later extended to the evaluation of proton 
strangeness form factors 




The calculation of Refs. [I^ll3|,|lj| did not remove the 
center-of-mass (CM) motion of the quark clusters. That 
could conceivably affect the estimate of the probability of 
the ss configuration from the measured strangeness mag- 
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netic moment considerably, as proposed in Ref. J12i |. Ac- 
cordingly, we set forth to reexamine the problem with the 
removal of the CM motion of the clusters. In the process, 
we obtained results which differ substantially from those 
presented in Refs. [12, llS uM where the CM motion was 
not removed. 

The configurations of the uudss component in the pro- 
ton considered in Refs. [12|, [ij, \l4\ are all of (4,1) clus- 



transition current to G^ 



tering type in that either four quarks uuds would be in P 
state with s in S* state (configuration I) or uuds in S state 
while s in P state (configuration II) , respectively, in order 
to ensure uudss as a whole has positive parity as proton. 
The symmetry of the spatial state of the four-quark system 
in configurations I and II would then be of [31] x and com- 
pletely symmetric [4]x, respectively. Within the harmonic 
oscillator constituent quark model, these two configura- 
tions are degenerate and of the lowest energy. However, 
the degeneracy is lifted by the color hyperfine quark-quark 
interaction as shown in Ref. [l^. After the splitting, the 
states of the lowest energy in configurations I and II for 
uuds cluster would have the flavor and spin state symme- 
try of [4]j;^s[22]_f[22]s and [31]fs[211]f[22]5, respectively 
12| . We will focus only on these two states of the lowest 
energy in this study. 

The calculation of the strangeness form factors G|;, G^^, 
and G\ within a harmonic oscillator constituent quark 
model is straightforward. However, the evaluation of the 
contribution of the transition current is involved if the CM 
motion is to be removed. 

We first compare our results for the form factors with 
those obtained in Ref. 13| where the CM motion is not 



removed. In order to facilitate the comparison, we will 
follow the notations and conventions of Refs. [l3|, |l4| as 
much as possible, unless otherwise specified. 

There are already some differences between our results 
and those given in Ref. [l3[ for the simple case of diago- 
nal matrix elements of the vector and axial vector current 
operators. Besides the factor of -^/l + q^ /Ami introduced 
in the denominators for G% and G\[ to account for rel- 
ativistic effects, we have not been able to reproduce the 
factor of (1 — (7^/18w|) (we use ws to denote the harmonic 
oscillator parameter associated with the uudss component 
so our a;5 is lo used in [l3|) which appears in the expres- 
sion for G\j. However, those differences are numerically 
insignificant in the low-q^ region where the nonrelativistic 
model calculations are, at best, expected to be valid. 

For the contribution of transition current, i.e., the non- 
diagonal matrix elements of the currents between 3g and 
Sg states, there are three discrepancies between our results 
and those of Ref. [13[. First, there is a sign difference 
for G|;, as well as that our result for G\j is -s/S times 
larger than that presented in Ref. [13|, which arises from 



a simple evaluation of the color matrix element |15| . The 
last and the most serious difference lies within G\ in which 
we obtain the following expression for the contribution of 
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where the superscript ND specifies that it is a non-diagonal 
matrix element. The parameters uj^, Puud and Ws, Pgs de- 
note the usual oscillator parameters and probabilities, re- 
spectively, of the uud and uudss configurations in the pro- 
ton, and G35 = (2w3a;5/(aj| +w|))^/^, while 5 denotes the 
relative phase between the uud and uudss components of 
the wave functions in the proton. Eq. ([T]) approaches con- 
stant at q^ = 0, while the corresponding result of Ref. 
13[ contains a factor of q^ and hence vanishes at q^ = 0. 
Since, as we see later, the contributions of non-diagonal 
matrix elements would dominate over the diagonal ones 
for all reasonable choices of wa and ws, our result of Eq. 
((T|) would lead to a consequence that G%j and G\ would 
be of the same sign, irrespective of the choice of phase 5, 
at low-q^ region which is in contradiction with the existing 
experimental data. It would then exclude the possibility 
that configuration I with four quarks uuds in P state and 
s in S' state could be the dominant configuration for the 
uudss, as concluded in Ref. [13[- We have carried out a 
direct numerical six-dimensional integration to verify that 
results agree with the analytical expression of Eq. ([IJ nu- 
merically. 

When the center-of-mass motion of the five-quark (5g) 
cluster is removed, we obtain the following results, in the 
case of configuration I, for the contributions of the diagonal 
matrix elements of the current to the proton strangeness 
form factors. 
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where the superscript D indicates that they are the di- 
agonal matrix elements. Without the removal of the CM 
motion of the 5(j cluster, the Gaussian factor e~ 
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Eqs. dHH) would become e"« /""= . This means that the 
removal of CM motion causes the form factors to decrease 
more slowly. Our results for the transition (3g — 5q) and 
(5(7 — 3q) contributions to the strangeness form factors, 
when the CM of 3q and 5q clusters are removed read as 



G'f'iq') = ^C.'/3^d^^L.-V/15^i 
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One observes that the exponents of the Gaussians in Eqs. 
([5]|31) and Eqs. dSJ?!) are different. This is expected since 
the center of masses of the 3g and 5q clusters are also dif- 
ferent. It is interesting to see that the transition current 
contributions drop faster than the diagonal ones. Further- 
more, it is seen that G^'^ and G^ are of the same sign 
as is the case before the removal of the CM motion, inde- 
pendent of the relative phase between the wavefunctions 
of 3g and 5g components. Consequently, as long as the 
transition current contributions dominate over the direct 
terms, then the configuration with s in 5 state cannot be 
the dominant configuration for uudss component. 

We next study configuration II which is degenerate 
with the lowest energy configuration I before being lifted 
by the color hyperfine quark-quark interaction. In config- 
uration II, uuds cluster is in S state while s is in P state. 
Only the results with the removal of CM motion will be 
presented. The direct current gives rise to the following 
contributions, 



G'f{q') - ll^"''^'"'^- 



We see that the direct current contributions to G\j and 
G\, as given in Eqs. (| MTU)) . are both negative at g^ = 
which contradicts the experiments. However, the tran- 
sition current contributions to G\,i and G\, as given in 

(■5) Eqs. (|12m3p are of opposite sign and in agreement with 
the data, independent of the sign of 5. Since the transi- 
tion current contributions dominate over the direct current 
contributions in the model considered here, it is hence of 

(6) interest to see whether wc could fit the experimental data 
of the proton strangeness form factors with configuration 
II or some linear combinations of configurations I and II. 
We take the proton and quark masses to be 0.938 and 

(.' / 0.313 GeV, respectively. The oscillator parameter for the 
3(7 core is fixed to be ^3 = 0.246 GeV. We then vary the 
oscillator parameter lo^ and the probability P^j of the 5g 
component uudss, to fit the experimental data G|; -I- rjG^j 
[8|, which are more directly measured in the experiments 
and G^ as extracted in Ref. [16[ where its sign in low-g^ 
region is well determined. Both signs of (5 = ±1 are tried 
and the best results arc then determined. 

Our best fits to the experimental data G% and G|; -|- 
TiG%,j, within configuration II, are shown in Fig. 1 as solid 
curves, with uj^ ~ 0.469 GeV, Psg = 0.025%, and relative 
phase Sp = +1, where subscript P refers to the fact that 
the configuration under consideration has s in P state. 
The data denoted with open circles are from Ref. [8[ and 
the crosses represent the corresponding values after the 
two-boson exchange effects are corrected [i7| . Solid trian- 
gles are the data from HAPPEx p and open boxes are 
from Ref. [l6|. The ensuing results for G|; and G^j are 
shown in Fig. 2. It is seen that the agreement with the 
data are in general quite good except for G|; + rjGfj and 
G|j at small values of g^, where there are large experimen- 
tal uncertainties. 

We have also explored the possibility of mixing config- 

\") urations II and 1, namely. 
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The results for the transition current matrix elements be- 
tween {iq — bq and bq — 3q) read as follows. 
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where \5q; a > and Sa denote the uudss states with s in 
cither 5 or P states and its relative phase with the three- 
quark state |3g >, respectively, to see whether a better 
description of the data can be obtained. It turns out that 
some improvements can be achieved only for G\ at low- 
q^ region with a small mixing probability of 6|. = 8% for 
configuration I, relatives phases Sp = l,Ss = —1, and a 
combined probability of Pss = ^5 = 0.058% (called ad- 
mixture A), as shown by the dashed curves in Figs. [Ijand 

m 

It is seen that we could fit the data reasonably well 
with a rather small probability of uudss component, e.g., 
Pss = 0.025% ~ 0.058% with either configuration II alone 
or a mixture of configurations II and I. It is in sharp con- 
trast to the values of Pss = 10 ~ 15% required in Ref. [13[ 
in order to fit G^^. Such a great reduction in Pss needed to 
reproduce the experimental data on the strangeness form 



Table 1: The strangeness magnetic dipole moment /Ib (in unit of /j]v), strangeness contribution to proton spin AS, and the radius of the uudss 
component, obtained in our calculation, with s in pure P state, S and P states admixture A, and B , as compared with the experiments. 





s in P state 
(0.025%) 


s in S and P admixture 
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Total 
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0.5 fm 
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0.096 
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Figure 1: Our predictions for form factors G| 



The 



results obtained with configurations with 5 in pure P state, S and P 
states admixture A, and B are denoted, respectively, by full, dashed, 
and dotted lines. Experimental data from Refs. Q, 0, llal are de- 
noted by open circles (GO), solid triangles (HAPPEx), and open 
boxes (Pate et al.), respectively, while the crosses are the corrected 
values of the GO data by taking into account the two-boson exchange 
mechanism [ITII. 



factors arises from several sources. These include the cor- 
rection factor v3 in the evaluation of color matrix ele- 
ments, change of the configuration for s from S* to P state, 
removal of the CM motion of the clusters, and the use of 
different model parameters. Each of them enhances the 
transition current matrix elements by ~^ 50 — 120%. It 
is interesting to note that our set of harmonic oscillator 
model parameters would give rise to a size of the uudss 
to be about 0.5 fm, which is quite close to that estimated 



by Rcf. [20| using a proton-core-^ picture for five-quark 
system with a scaling factor s = 1.5. 

The corresponding results for strangeness magnetic mo- 
ment fXs, strangeness contribution to proton spin AS' = 
6*^(0), and the size of uudss component r^q, with s in 
P state or admixture of S and P states, are given in Ta- 
ble 1 and compared with the experiments. The agreement 
between experiments and numbers obtained with s in ad- 
mixture of S and P states seems reasonable. 
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Figure 2: Our results for G"^, and G\j. Notations are the same as 
in Fig. [T] Experimental data are from Refs. [3, 0, 0, |a and denoted 
with solid diamonds (SAMPLE), solid triangles (HAPPEx), solid 
boxes (A4), and open circles (GO), respectively. 



Recently, Chang and Peng [2l| generalize the approach 
of Brodsky, Hoyer, Peterson, and Sakai (BHPS) [22| for 
the intrinsic charm quark distribution in the nucleons to 



the light-quark sector involving intrinsic u, d, and s sea 
quarks to analyze the existing d{x) — u{x), s{x) + s(x), and 
u{x) + d{x) — s{x) — s{x) data and conclude that probabil- 
ity Pss for five-quark configuration uudss to lie between 
0.024-0.029. We explore the consequence of such a result 
to our model calculation by fixing P^s = 2.4% and vary uj^ 
to fit the data. The resultant fit we obtain with tu^ = 0.108 
GeV, which corresponds to a large size of the five-quark 
system with r^q = 2.16 fm, and a small admixture of S 
state with a probability of about 15% (called admixture 
B), are shown in Figs. [1] and [5] by dotted lines. The most 
interesting feature of this fit is the appearance of a bump 
in C'^+rjC'^.j in the very low-q^ region with q^ < 0.1 GeV^, 
which seems to be hinted by the GO data but hampered by 
large experimental error bars and fluctuations. It would be 
worthwhile to carry out experiments in such a low-q^ re- 
gion if further theoretical study would support this behav- 
ior. All form factors vanish rapidly beyond q^ > 0.2 GcV^ 
because of the small value of a;5 and Gaussian nature of 
the harmonic oscillator wavcfunctions. Wc note that the 
predicted values of fig, A5, and rsg, as presented in Table 
1, agree with the data within experimental errors. How- 
ever, the experimental value of 0.37 ± 0.79 was obtained 
in Ref. [18[ in a fit to the form factors including the ncxt- 
to-leading-order terms of the Q^- dependence. In another 
fit with only leading order in Q^ which has almost same 
reduced x^i a value oi fig = 0.12 ±0.55 ±0.07 was obtained 
which would disfavor configuration B. 

In summary, we have reinvestigated, within a nonrela- 
tivistic constituent quark model, the question of whether 
a five-quark component with configuration of (4,1) clus- 
tering, as previously considered by Riska and Zou [12l Il3j 
can account for the data of the proton strangeness form 
factors. Two configurations of the lowest energies, both 
consist of four quarks in colored state and one antiquark 
are considered. They possess spatial-fiavor-spin symmetry 
[31]x-[^]fs[22]f[22]s and [4]^- [31]f5[211]f[22]s , with 
antiquark s in the S (configuration I) and P states (config- 
uration II), respectively. They are degenerate before be- 
ing splittcd by the color hyperfine quark-quark interaction 
with configuration II of higher energy. 

Wc have not been able to reproduce the results of Ref. 

13l | to substantiate their claim that configuration I is to be 

preferred as the dominant configuration of a pos sible five- 



quark component. The claim made in Ref. [13| is based 
on a particular choice of the relative phase S between uud 
and uudss components of the proton wave functions such 
that G|j could be positive at small q^ region. However, 
the corrected expression of G^ we obtain for configura- 
tion I would lead to a G\ to be of the same sign as G\j 
in the low-q^ region which clearly contradicts all existing 
data. Wc also find that the expression for the transition 
current contribution to G|; obtained in Ref. [l3| carries 
an erroneous minus sign. 

We then proceed to study configuration II, where s is 
in P state, and make an effort to remove the CM motions 
of the clusters which was not done in Ref. (iSi . We demon- 



strate that it is possible to give a satisfactory description 
of the existing data on the proton strangeness form factors 
G|;, G^j, G\, and the linear combination of G|; ± 'qG\j 
which is more directly extracted from the parity-violating 
asymmetry Apy measured in elastic electron-proton scat- 
tering with a very small probability of Pss = 0.025% for 
the uudss components. The agreement with G\ data can 
be improved in the low-g^ region by considering an admix- 
ture of configurations I and II with a total uudss probabil- 
ity Pss increased to 0.058% with configuration I accounts 
for 8% of the total. Wc further find that without removing 
CM motion, Pss would be overestimated by about a factor 
of four in the case when transition current dominates. 

We have also explored the consequence of a recent 
claim [2l| , reached from analyzing existing data on d{x) — 
u{x), s{x) + s{x), and u{x) + d{x) — s(x) — s{x), that Pss 
lies between 2.4— 2.9%. A small bump in both G% + r]G%j 
and G|; in the region of q^ < 0.1 GcV^ for an admixture 
of configuration I and II. 

It is tempting to conclude, from the results presented 
above, that the current strangeness form factors experi- 
ments seem to indicate that the dominant configuration 
in which uudss arrange themselves in configuration II, in 
which uuds is in a completely symmetric spatial state [4]x 
with fiavor and spin state symmetry of [31]i^5[211]i?[22]s. 
However, there are a few caveats here. First is that the 
agreement between our results and the existing data is not 
perfect, to say the least. For example, recent data from 
A4 at q^ = 0.22 GcV^ gives a negative value of G\^ = 
—0.14 ±0.11 ±0.11. One might also ask whether a nonrel- 
ativistic constituent quark model is quantitatively reliable 
in evaluating the contributions of transition current which 
is found to be dominant in our calculation but is of a rel- 
ativistic effect in nature. More study about this issue is 
clearly needed. Lastly, it is known that configurations of 
(4, 1) clustering and (3, 2) clustering, like the meson cloud 
configuration, are not orthogonal. Then whether (4, 1) or 
(3, 2) configuration would be favored might be just a mat- 
ter of choosing a different basis if it turns out that more 
precise data will require linear combination of several (4, 1) 
clustering configurations. 
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